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Abstract 

We present a new computation of S'-wave color-singlet nonrelativistic QCD matrix elements for 
the J ftp and the rjc- We compute the matrix elements of leading order in the heavy-quark velocity 
V and the matrix elements of relative order v"^. Our computation is based on the electromagnetic 
decay rates of the J/ip and the r]c and on a potential model that employs the Cornell potential. We 
include relativistic corrections to the electromagnetic decay rates, resumming a class of corrections 
to all orders in v, and find that they significantly increase the values of the matrix elements of 
leading order in v. This increase could have important implications for theoretical predictions for 
a number of quarkonium decay and production processes. The values that we find for the matrix 
elements of relative order v'^ are somewhat smaller than the values that one obtains from estimates 
that are based on the velocity-scaling rules of nonrelativistic QCD. 
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I. INTRODUCTION 



In the nonrelativistic quantum chromodynamics (NRQCD) factorization formalism [1], 
heavy-quarkonium decay and production rates are expressed as sums of short-distance coef- 
ficients times NRQCD operator matrix elements. The matrix elements in these sums scale 
as powers of v, the typical heavy-quark (or antiquark) velocity in the quarkonium rest frame. 
Hence, the sum in the NRQCD factorization expression can be thought of as an expansion 
in powers of v. The term that is proportional to the matrix element of leading order in v 
often gives the dominant contribution in decay and production processes. The leading-order 
matrix element involves the production or annihilation of a heavy quark- antiquark (QQ) 
pair in a color-singlet state. The term that is proportional to the matrix element of relative 
order v"^ gives the first relativistic correction. This order-v^ matrix element also involves the 
production or annihilation of a heavy QQ pair in a color-singlet state. 

In the vacuum-saturation approximation [1] for decay matrix elements, one keeps only the 
vacuum intermediate state, while, in the vacuum-saturation approximation for production 
matrix elements, one keeps only the heavy QQ intermediate state. The vacuum-saturation 
approximation is valid up to corrections of relative order v"^ [1]. In this approximation, the 
color-singlet decay matrix elements are equal to color-singlet production matrix elements. 
These vacuum-saturation matrix elements arc also the relevant ones for purely electromag- 
netic decay and production processes and for exclusive decay and production processes. The 
vacuum-saturation matrix element at leading order in v is proportional to the square of the 
quarkonium wave function at the origin. In this paper, we compute the vacuum-saturation 
matrix elements of leading order in v and of relative order v'^ for the J/ip and rjc states. 

The analysis of these matrix elements for the J/i/j state differs in several respects from 
a previous one involving some of the authors [2]. In that analysis, the matrix element at 
leading order in v was obtained by comparing the theoretical expression for the decay width 
r[J/ip e+e"] with the experimental measurement. The theoretical expression that was 
used in that analysis included the order-cts correction, but not the relativistic corrections. In 
the present paper, we include those relativistic corrections. The matrix element of relative 
order is determined from a potential-model calculation [2] that uses the leading-order 
matrix element as an input. The relativistic corrections to T[J/iIj e+e^] in turn depend 
upon that order-v^ matrix element. Hence, the leading-order and order-v^ matrix elements 
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are related through a coupled pair of nonlinear equations, which we solve numerically. 

We obtain values for the rjc matrix elements in two different ways and average the results. 
First, we obtain a set of values by making use of the comparison between theory and experi- 
ment for the width T[r]c — > 77]. This comparison gives one nonlinear equation for the matrix 
elements. As in the J/ ip case, we make use of a potential-model calculation of the order-v^ 
matrix element to obtain a second nonlinear equation, and we solve the coupled nonlinear 
equations numerically to obtain a set of values for the rjc matrix elements. We obtain a 
second set of values by making use of the fact that, because of the approximate heavy-quark 
spin symmetry of NRQCD [1], the rjc and J/ip matrix elements are equal, up to corrections 
of relative order v^. Wc define this second set of vahics for the ric matrix elements simply 
by taking the values that we obtain for the J/ip matrix elements and appending additional 
error bars that take into account the order-v^ corrections to the heavy-quark-spin-symmetry 
relation. 

Because the two sets of values for the rjc matrix elements that we obtain in this way have 
input parameters (such as the heavy-quark mass and the string tension) in common, the 
uncertainties in these matrix elements are highly correlated between sets and between matrix 
elements within a set. Therefore, we carry out a covariance-matrix analysis to compute the 
average. The J/ip and r}c matrix elements are also highly correlated. Such correlations could 
be important in applications of our results to calculations involving both the J/ip and the 
Tjc and/or order-v^ corrections. Therefore, we present tables showing the variations of each 
of the matrix elements with respect to the various sources of uncertainty and also give the 
covariance matrix that corresponds to these variations. 

A further refinement that we include in this work is to resum a class of relativistic 
corrections [2] to r[J/ijj e^e~] and T[ric 77]. First, we consider all corrections that 
arise from matrix elements involving only color-singlet QQ Fock states. By making use of a 
generalization of the Gremm-Kapustin relation [2, 3], we can determine all of these matrix 
elements, up to corrections of relative order from the leading-order and order-i>^ matrix 
elements. The simple expressions that result can easily be summed to all orders in v. This 
resummation is equivalent to retaining all of the relativistic corrections that are contained 
in a potential-model QQ wave function, up to the ultraviolet cutoff of the NRQCD matrix 
elements. 

Because the expressions for the matrix elements of order v"^ and higher are accurate only 
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up to corrections of relative order v^, the uncertainty in the resummed expression is of order 
v'^ relative to the leading-order expression. That is, the nominal accuracy in v is no higher 
than that of a fixed-order calculation through relative order v^. However, if the relativistic 
corrections to a given process that arise from the QQ Fock-state wave function have partic- 
ularly large coefficients in the v expansion, then the use of the resummed expression may 
improve the numerical accuracy. Furthermore, the resummation may give an indication of 
the rate of convergence of the v expansion. In any case, it is generally desirable to include 
in a calculation a well-defined, if incomplete, set of contributions whenever possible. 

The remainder of this paper is organized as follows. In Sec. II, we review the definitions 
of the iS-wave NRQCD matrix elements at the leading and higher orders in v, and we give 
the relations of these matrix elements to the quarkonium wave functions. We also introduce 
the generalized Gremm-Kapustin relation for an 5'-wave quarkonium state, which expresses 
matrix elements of higher order in v in terms of the matrix element of leading order in v and 
the binding energy. The generalized Gremm-Kapustin relation allows us to resum a class of 
relativistic corrections to quarkonium decay to all orders in v. In Sec. Ill, we present the 
resummed formulas for the electromagnetic decay widths of the J/if) and the rjc- Sec. IV 
contains a description of the potential-model method that we use to compute the binding 
energy of the 5'-wave states and, through the generalized Gremm-Kapustin relation, the 
NRQCD matrix elements of higher order in v. In Sec. V, we compute the numerical values 
of the NRQCD matrix elements for the J/t/j and the rjc- We compare our results for the 
matrix elements with those from previous determinations in Sec. VI. Finally, we summarize 
our results in Sec. VII. 

II. NRQCD MATRIX ELEMENTS 

A. Decay and production matrix elements 

In the cases of the inclusive decays of spin-singlet and spin-triplet 5'-wave quarkonium 
states, such as the rjc and the J/ip, the matrix elements at the leading power in v are 

{0,CSo))h = {HCSo)\^hx^mCSo)), (la) 

m's,))H = (//e^oivVx • xVV'i/zf (lb) 
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where if is a quarkonium state, '^^'^^Sj is the standard spectroscopic notation for a state with 
spin angular momentum s, orbital angular momentum zero, and total angular momentum 
J, is a two-component Pauli spinor that annihilates a heavy quark, x is a two-component 
Pauli spinor that creates a heavy antiquark, and is a Pauli matrix. The subscript 1 on 
an NRQCD operator O indicates that it is a color-singlet operator. 

Similarly, in the case of the inclusive production of spin-singlet and spin-triplet S'-wave 
quarkonium states, the matrix elements at the leading power in v are 

= (0|xV I E \HCSo)+X){HCSo)+X\U^x\0), (2a) 
\x, pol. / 

m's,)H) = (o|xV> ( \HCS,) + X){H{'S,) + X\ J ^V^xlO), (2b) 

\x, pol. / 

where the sum is over the light degrees of freedom X and the 2 J+1 quarkonium polarizations. 

In the vacuum-saturation approximation [1], which is valid up to corrections of relative 

order the decay matrix elements in Eq. (1) and 1/(2 J -|- 1) times the production matrix 

elements in Eq. (2) both reduce to 

(Oi(%))^s = |(0|xVl^(%))r, (3a) 

= mx^a-e'^mCSiW, (3b) 

in the spin-singlet and spin-triplet cases, respectively. In Eq. (3b), the quarkonium po- 
larization vector is denoted by e, and there is no sum over the polarization states of the 
quarkonium. In the cases of purely electromagnetic decay or production or exclusive decay 
or production, the matrix elements in Eq. (3) are the relevant ones at leading order in v. 

The first relativistic corrections to inclusive S-wave decay and production involve oper- 
ators that are analogous to those in Eqs. (1) and (2), but that contain a factor of (— 1£))^ 
between either tp^ and x or x^ and jp. Here, £> is the spatial part of the covariant derivative 
acting to the left and right anti-symmetrically: x^-^V' = xH^i') ~ i^xVi'- These operators 
are of order v'^ relative to those in Eqs. (1) and (2). The corresponding matrix elements 
reduce in the vacuum-saturation approximation to 

{riCSo))l' = Re [{HCSoM^x\0){0\x\-|Dr^|;\HCSo))], (4a) 
{rii'S,))l' = Re [{Hi'S^M^a- • ex|0)(0|xV • e*(-iB)V|//(=^50)] ■ (4b) 



In the cases of purely electromagnetic production or decay or exclusive production or decay, 
the matrix elements in Eq. (4) are the relevant ones. 

Corrections of still higher orders in involve, among other matrix elements, those in 
which higher powers of (— |Z))^ appear. It is convenient to define ratios of these matrix 
elements to the matrix elements of lowest order in v : 
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where q is half the relative three-momentum of the Q and Q in the quarkonium rest frame. 

In this paper, we compute the quantities {Oi{^ SQ))^f and {Oii^ Si))^jJ^, which are given 
by Eq. (3), and the quantities {q'^^)r]c (9^'").//V'' which are given by Eq. (5). As we shall 
see, the higher-order ratios in Eq. (5) can be related to the lowest-order ones by making use 
of a generalization of the Gremm-Kapustin relation [2, 3]. 

As is discussed in Ref. [2], the higher-order matrix elements in Eq. (5) contain power 
ultraviolet divergences and require regularization. In this paper, we regulate these power 
divergences dimensionally at the one-loop level. One-loop dimensional regularization of 
the matrix elements is appropriate for use in conjunction with one-loop calculations of the 
short-distance coefficients. 



B. Relations of NRQCD matrix elements to quarkonium wave functions 

In the rest frame of an jS'-wave heavy quarkonium H in a spin-singlet C'Sq) or spin-triplet 
{^Si) state, one can express the wave function at the origin of the leading QQ Fock state in 
terms of the color-singlet NRQCD matrix elements [1]: 

M^sd^) = J -0^sM^s,){q) = ^(0|xV|/^(^^o)), (6a) 

6^^^(350(0) = 6 J ^,M^s,){q) = -^m^rrmeSi)). (6b) 

iPh{q) is the momentum-space wave function for the leading Q{q)Q{—q) Fock state of the 
quarkonium. The wave function is, of course, gauge dependent. Throughout this paper, we 
work in the Coulomb gauge. The normahzation factor l/y/2Nc accounts for the traces in 
the SU(2)-spin and SU(3)-color spaces. 
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Relativistic corrections to the production and decay rates for a heavy quarkonium involve 
matrix elements that are related to derivatives of the wave function at the origin: 

<i5o)(0) = / ^,<l''M^s,M) = ^(0|x^(-iV)^>|i/(^5o)), (7a) 



e 



Usually, these operator matrix elements are written in terms of the covariant derivative D 
(Ref. [1]), as in Eqs. (4) and (5), rather than V. However, in the Coulomb gauge, the 
difference between the D and V is suppressed as (Ref. [1]). We emphasize again that 
the derivatives of the wave function at the origin, as defined in Eq. (7), are ultraviolet- 
divergent quantities, which must be regulated. We also note that different from 
the expectation value of g^^: 

7^/(03 q''rHmH{q). (8) 

Comparing Eq. (7) with Eq. (5), we see that 

One can also define matrix elements of powers of the heavy-quark velocity in terms of matrix 
elements of powers of the heavy-quark momentum: 

{v'')h = (10) 
where mq is the heavy-quark mass. 



C. The generalized Gremm-Kapustin relation 

Prom the effective field theory known as potential NRQCD (pNRQCD) [4] , it follows that 
one can compute the wave functions at the origin and derivatives of wave functions at the 
origin in Eq. (9), up to errors of relative order f^, from the Schrodingcr wave function for a 
heavy QQ pair interacting through the leading (static) spin-independent QQ potential. It 
was shown in Ref. [2] that, in the case of a spin- independent potential and for dimensionally 
regulated matrix elements, the ratios in Eq. (9) are related through the generalized Gremm- 
Kapustin relation: 

{q'^)H = {me^sr[l + 0{v')l (11) 
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where CnS is the binding energy of the QQ pair in the quarkonium state H with principal 
quantum number n and orbital angular momentum and m is the heavy-quark mass in 
the effective theory pNRQCD. The relation (11) follows from the equations of motion of 
the QQ pair and from dimensional rcgularization of the matrix elements at the one-loop 
level, provided that the potential is parametrized as a sum of constants times powers of 
the QQ separation. The Cornell potential [5], which we will employ later in this paper, is 
parametrized in this way. We note that Eq. (11) implies that 

{q'lH = {q'Tu. (12) 

up to corrections of relative order v'^ . 

We will use Eq. (11) to determine the quantities {q'^^)nc {q^^)j/ti>- order to evaluate 
the ground-state binding energy 615, we will make use of a potential model that is based on 
the Cornell potential. 

III. FORMULAS FOR ELECTROMAGNETIC DECAYS OF S-WAVE HEAVY 
QUARKONIA 

In this section, we present the NRQCD factorization expressions for the electromagnetic 
decay widths r[i7(^5'i) e+e"] and r[i7(^S'o) 77]. In subsequent parts of this paper, 
we will apply these formulas to the decays J/%1) ^ e^e~ and r]c — > 77. 

A. V[H{^S^) e+e-] 

The NRQCD factorization formula for the amplitude for the decay H(J'Si) — » e+e~ is 

A[H{'S^) ^ e+e-] = v^5^rf„('5i)(0|On|i/(^5i)), (13) 

n 

where ttih is the quarkonium mass, the dn{^Si) are short-distance coefficients, and the On 
are NRQCD operators. The prcfactor y/2mH compensates for the fact that the hadronic 
NRQCD operator matrix elements conventionally have nonrelativistic normalization, while 
we choose the amplitude on the left side of Eq. (13) to have relativistic normalization. 

Now we approximate the formula (13) by retaining only those operator matrix elements 
that connect the vacuum to the color-singlet, QQ Fock state of the quarkonium H. Then, 
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we have 

n 

where the short-distance coefficients c^(^5'i) are a subset of the short- distance coefficients dn 
in Eq. (13). We will clarify the meaning of the approximation that we have taken to arrive 
at Eq. (14) below. 

Because the c^(^5'i) are insensitive to the long-distance nature of the hadronic state, we 
can calculate them by replacing the initial hadronic state \^2mH\H{^Si)) in Eq. (14) with 
a perturbative spin-triplet S'-wave QQ state: 

Am^eS,) ^ e+e-] = 5^<(^5i)(0|x^(-|B)2"<7V|Qgi(^5i)). (15) 

n 

The factor ^2mu is absent in Eq. (15) because we use the same (relativistic) normalization 
for the QQ state on both sides of Eq. (15). 

In the rest frame of the quarkonium, the perturbative amplitude on the left side of Eq. (15) 
at order q;° is [6] 

^ (16) 



2 / 

^[ggi('^i) ^ e+e-] = .JW,2E{qf-^ L ■ e[l - 



3E{q)[E{q)+mQ]J' 

where e is the electromagnetic coupling constant, eg is the electric charge of the heavy quark, 
e is the polarization vector for the spin-triplet state, L is the leptonic current, and E{q) = 



itlq + q^. In the expression (16), wc have neglected the electron mass in comparison 
with the quarkonium mass. The factor l/rn\j arises from the photon propagator. The 
perturbative matrix elements on the right side of Eq. (15) are given by 

(o|x^(-|B)2vv|ggi(35i)) = ^,2E{q) ^v. (i?) 

The factor 2E{q) arises from the relativistic normalization of the QQ state. By comparing 
Eqs. (15) and (16), one can read off the short-distance coefficients c\^{^Si)\ 



j(3o) ^ ^ r (A_Y (l ^ 



(18) 

q2=0 



Substituting the c\^{^Si) in Eq. (18) into Eq. (14), using Eq. (12), and including the 
order-Qis correction to the amplitude [7, 8] , we obtain 

A[H{^S,) e+e-] = ^/^^^t^V \l - /((q^)^^/m^) - 2Cf-] {Q\x^a^^P\H{^S^)), (19) 



m 



H 



TT J 



where f{x) is defined by 

3{1 + X + y/l + X) 

Now we can clarify tlie meaning of the approximation that was taken to arrive at Eq. (14). 
Suppose that we specialize to the Coulomb gauge. Then, we can drop the gauge fields in 
covariant derivatives in the matrix elements in Eq. (14), making errors of relative order 
v^. The matrix elements are then proportional to derivatives of the Coulomb-gauge color- 
singlet QQ quarkonium wave function at the origin [1]. (See Sec. II B.) That is, they are 
proportional to the moments of the momentum-space wave function with respect to the 
wave-function momentum (the relative momentum of the Q and Q). From Eq. (18), we 
see that the short-distance coefficients cl^{^Si), when contracted into e\ are the coefficients 
of the Taylor expansion of A[QQi{^Si) — > e'^e~]/[y/2Nc2E{q)] with respect to the wave- 
function momentum. Hence, Eq. (14) has the interpretation of the convolution of the short- 
distance amplitude with the momentum-space quarkonium wave function, where the short- 
distance coefficients have been Taylor expanded with respect to the wave-function momenta. 
Therefore, we sec that the approximate NRQCD expansion in Eq. (14) includes all of the 
relativistic corrections that are contained in the color-singlet QQ quarkonium wave function, 
up to the ultraviolet cutoff of the NRQCD matrix elements.^ 

We note that, in the quarkonium rest frame, the square of the spatial part of the leptonic 
factor L, summed over lepton spins, is given by 

J2 L'L*^ = 2m| (S'^ - k^y^ , (21) 

spins 

where k — k/\k\ and k is the three-momentum of the e~ in the quarkonium rest frame. The 
temporal parts of Jospins L^L*'^ vanish in the quarkonium rest frame. 

We obtain the leptonic decay width of the spin-triplet S'-wave heavy quarkonium by 
taking the square of the amplitude (19), summing over lepton spins using Eq. (21), averaging 
over the H(^Si) polarization states, and multiplying by the two-body phase space and the 
normalization {2mH)~^- The result is [6, 9] 

T[HCS,) ^ e+e-] = [l - fi{q')H/ml) - 2Cj.^ ' {0,)h, (22) 



TT J 



^ We note that, in the case of dimensionaUy regulated NRQCD matrix elements, pure power ultraviolet di- 
vergences in the matrix elements are set to zero. Hence, the effects of purely power-divergent contributions 
are absent in the resummation. 
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where a — /{At^). In Eq. (22), the exphcit relativistic corrections are contained in the term 
— f {{q^ / mq) h) ■ In addition, there are imphcit relativistic corrections that are contained in 
the factors rriH- Strictly speaking, if one were to compute the decay amplitude completely 
with the framework of NRQCD, then ttih would be written as 2E{q) and expanded in powers 
of \q\/mQ to obtain additional relativistic corrections. (See, for example, Refs. [6, 10].) 
However, we note that the factor l/rrijj in Eq. (22) is clearly identifiable as arising from the 
photon propagator and the leptonic current, and, so, it is not necessary to treat that factor 
within the framework of NRQCD. We choose not to apply the nonrelativistic expansion of 
NRQCD to the factor l/mjj. That is, we apply NRQCD only to the heavy-quark factor in 
the amplitude. This choice reduces the theoretical uncertainties by making use of the fact 
that the quarkonium masses are known very precisely. 

The order-a^ corrections to r[i7(^5'i) —>■ e+e"] (Ref. [11, 12]) contain a strong dependence 
on the NRQCD factorization scale. If one were to include those corrections in the expression 
(22) and use it to determine then {Oi)h would also contain a strong dependence 

on the NRQCD factorization scale. If one were to make use of {Oi)h in calculating other 
quarkonium decay and production processes, then the factorization-scale dependence would 
cancel only if the short-distance coefficients for those processes were calculated through 
relative order a^. Generally, short-distance coefficients for quarkonium processes have not 
been calculated beyond relative order a^. For this reason, we have chosen to omit the 
order-a^ corrections to the leptonic width in Eq. (22). 

B. r[/f(i5o) ^ 77] 

Employing a method analogous to that which is given in Sec. Ill A, one can obtain the 
NRQCD factorization formula for the relativistic corrections to the two-photon decay of a 
spin-singlet 5'-wave quarkonium state H{^So). 

The NRQCD factorization formula for the amphtude for the decay H(^So) — > 77 is 

A[HCSo) ^ 77] = V2^J2'^r^CSomOr,\HCSo)), (23) 

n 

where the dn{^So) are short-distance coefficients. We approximate this expression by keeping 
only those matrix elements that connect the vacuum to the color-singlet QQ Fock state in 
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the quarkonium. Then, we have 



(24) 



As in the spin-triplet case, this modified NRQCD factorization formula retains all of the 
relativistic corrections that are contained in a potential-model QQ wave function, up to the 
ultraviolet cutoff of the NRQCD matrix elements. 

We can calculate the short-distance coefficients c„(^S'o) by replacing the initial hadronic 

state ^/2mH\H{^So)) in Eq. (24) with a perturbative spin-singlet S'-wave QQ state: 



A[QQiCSo) ^ 77] = 5^c„(i5o)(0|x^(-iB)^>|QQi(^5o)). 



(25) 



In the rest frame of the quarkonium, the perturbative amplitude on the left side of Eq. (25) 
at order a° is given by [6] 

MQQ.CS,) ^ = V^e^eJ m (26) 

\q\ J^[q) - \q\ 

where ki and are the momentum and the polarization of the i-th photon. The perturbative 
NRQCD matrix elements on the right side of Eq. (25) are given by 



(0|x^(-|r>)^"V'|QOi(%)) = A/2A^2£;(5)q2". 



(27) 



By comparing Eqs. (25) and (26), one can read off the short-distance coefficients c„(^5'o): 



CnCSo) = e^e^ 



fci • €1 X e*2 

\ki\ 



d 



(28) 



n\ \dqy 2E{q)\q\ E(q) - \q\ 

Substituting the c„(^5'o) in Eq. (28) into Eq. (24), using Eq. (12), and including the 
order-CKg correction to the amplitude [13-15], we obtain 

20 - 



A[H{ So) 77J = e Cr 



mq |fei| 
x(0|xV|^(%)), 



1 - 9 {{q^)H/ml) 



-Up — 



TT 



(29) 



where g{x) is defined by 

g{x) = 1 
= 1 



1 



2^a;(l + a;) 
1 

2^a;(l + a;) 



log 



log 



1 + 2^/x{lTx) + 2x 



(30) 
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Squaring the amplitude (29), summing over the photon polarizations, multiplying by the 
phase space and 1/2! for the two identical particles in the final state, and dividing by the 
normalization 2mH, we obtain the two- photon decay width of a spin-singlet 5"- wave heavy 
quarkonium: 

2 

{Oi)h. (31) 

In the formula (31), we have omitted the order-a^ corrections to the decay amplitude 
[16]. As we discussed in the case of the leptonic width of a spin-triplet 5"- wave quarkonium, 
the order-a^ corrections contain a dependence on the NRQCD factorization scale and can 
only be used consistently in conjunction with calculations of other quarkonium processes 
through order a^. 

IV. POTENTIAL MODEL 

As we have explained earlier, in order to compute the higher-order matrix elements that 
appear in Eq. (5), we need to compute the ground-state binding energy eis that appears 
in the generalized Gremm-Kapustin relation (11). In this section, we describe briefly the 
potential model that we use to compute eis- For details of the model, we refer the reader 
to Refs. [2, 5]. 

The model makes use of the Cornell potential [5], which parametrizes the QQ potential 
as a linear combination of the Coulomb and linear potentials: 

V{r)^--+ar, (32) 
r 

where k is a dimensionless model parameter for the Coulomb strength and a is the string 
tension, which is of mass dimension two. In the original formulation of the Cornell potential 
model [5], the strength of the linear potential was given in terms of a parameter a, where 

a^l/^/a. (33) 

By varying the parameters in the Cornell potential, one can obtain good fits to lattice 
measurements of the QQ static potential [17]. Therefore, we assume that the use of the 
Cornell parametrization of the QQ potential results in errors that are much less than the 
order-i>^ errors (about 30%) that are inherent in the leading-potential approximation to 
NRQCD. 
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mCSo) ^ 77] - ^ 



2\ / 2^ 20 -TT^ as 

~" -Up — 



l-g{{q')H/ml)- 
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The Schrodinger equation for the radial wave function Rneir) with the radial and orbital 
angular-momentum quantum numbers n and i is 

1 d f o d 



(34) 



where m is the quark mass and is the binding energy of the n£ state. We treat m as a 
parameter of the potential model and note that it is, in general, different from the heavy- 
quark mass niQ, which appears in the short-distance coefficients of NRQCD factorization 
formulas. As usual, for an 5"- wave state, the wave function is ■0ns('") — R-ns{''^)/V^- 
Introducing the scaled radius p and scaled coupling A [5], 



p — {amY^^ r, 

A = ^ 

((7/m2)V3' 

which are dimensionless, one can rewrite the radial equation (34) as [5] 

dp'^ " 



+ 1 A 

+ --P + Cn 
p^ p 



(35a) 
(35b) 



(36) 



where Uniip) Cnf. arc the dimensionless radial wave function and the dimensionless energy 
eigenvalue of the n£ state. The relation between Rni{r) and Une{p) is 

Une{p) 



Rntir) 



'am ■ 



where the wave functions are normalized according to 

poo poo 

/ \uM\Hp^ / |i?„,(r)|V2dr = l. 

JQ Jo 

The binding energy is related to the dimensionless eigenvalue as 



(37) 



(38) 



(39) 



Now let us specialize to the 5"- wave case. In order to compute CnS from Eq. (39), we 
must fix the model parameters a, m, and A and solve Eq. (36), with i = 0, for (nsW- Our 
strategy is to fix a from lattice measurements and to use the measured 1S-2S mass splitting 
and |V'n5(0)P, as determined from the electromagnetic decay widths, to solve for m and A. 
Using Eq. (39), we can express m in terms of the 1S-2S mass sphtting: 



m(A) 



C25(A) - Cig(A) 

77125 - mis 



(40) 
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For S-wsive states, the wave function at the origin ■0n5(O) = Rns{0) / can be expressed 
as [5] 

l^nsm' = £ / d\\^r,s{r)\'^ = ^ [1 + AF„,(A)] , (41) 
where F„5(A) is the expectation value of for the nS state: 

FnsiX)^ r% \Uns{p)?. (42) 

Jo P 

The first equahty in Eq. (41) can be obtained by multiplying the radial Schrodinger Equation 
(34) on the left by R^g{r) and integrating by parts. 

For purposes of computation of the NRQCD matrix elements, it is convenient to express 
those matrix elements in terms of the potential-model parameters. From Eqs. (39) and (41) 
and the generalized Gremm-Kapustin relation (11), we find that 

{0,)h = 2N,\m\' = [1 + XF,s{X)] , (43a) 

{q')H = m{\)e,s{\) = [am{\)]'/\,s{X): (43b) 

where m(A) is given in Eq. (40) and Fis{X) is given in Eq. (42). 

V. COMPUTATION OF THE NRQCD MATRIX ELEMENTS 

In this section, we determine the numerical values of the NRQCD matrix elements for the 
J/ip and the r]c- In this and subsequent discussions, we drop the superscript VS on {Oi)^ 
because other sources of uncertainty, which we will describe, are much larger than the error 
in the vacuum-saturation approximation. 

A. Method of computation 

Were it not for the relativistic corrections in the decay widths in Eqs. (22) and (31), 
we could simply solve those equations for {Oi)h- Then we could use the value for {Oi)h 
that we would obtain to solve Eq. (43a) for A and use that value of A to solve Eq. (43b) 
for (qf^)ij . Because the relativistic corrections in Eqs. (22) and (31) couple those equations 
weakly to Eq. (43b), we must carry out the more difficult task of solving Eq. (22) or Eq. (31) 
simultaneously with Eqs. (43a) and (43b). 
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First, we express r[J/ip — > e"^e~] and r[77c — > 77] in terms of the potential-model pa- 
rameters by substituting Eq. (43) into Eq. (22) and Eq. (31), respectively. We equate those 
expressions to the experimental values of the electromagnetic widths [18]: 

^[J/^jJ e+e"] = 5.55 ± 0.14 ± 0.02 keV, (44a) 
r[r7c ^ 77] = 7.2 ± 0.7 ± 2.0 keV. (44b) 

Then, we solve the resulting equations numerically for the model parameter A. In computing 
the solution, we express the eigenvalues Cis{^) and C2s{^) and the expectation value Fis{X) 
[Eq. (42)] as functions of A by fitting interpolating polynomials to computations of the 
eigenvalues and expectation value at fixed values of A. Once we have obtained a value for 
A, we substitute it into Eq. (43) to obtain values for the NRQCD matrix elements. 

In carrying out the numerical computation, wc need values for the charm-quark mass rric, 
the string tension a, and the 1S-2S mass splitting. In order to maintain consistency with 
the calculations of the electromagnetic decay widths of the J/t/j and the rjc at NLO in cts, 
we take rric to be the pole mass. The specific numerical value that we use is^ 

rric = 1.4 ± 0.2 GeV. (45) 

We fix the string tension a by making use of lattice measurements. Prom Ref . [19] , we find 
that aal = 0.0114(2) at a lattice coupling j3 = 6.5, where is the lattice spacing. Lattice 
calculations of the hadron spectrum at = 6.5 yield values for I/ol of 3.962(127) GeV 
(Refs. [20, 21]) and 3.811(59) GeV (Refs. [20, 22]). These result in values for the string 
tension of cr = 0.1790±0.0119 GeV^ and a = 0.1656±0.0059 GeV^, respectively. Combining 
these two values, we obtain 

(T = 0.1682 ± 0.0053 GeV^ (46) 
For the 1S-2S mass splitting, we take the mass difference between the J/ip and ip{2S) [18]: 

m2s - mis = 589.177 ± 0.036 MeV. (47) 

We use mj/^ — 3.096916 GeV and m^^ = 2.9798 GeV [18]. We also need values for ag. In 
the case of J/^ — > e"^e~, we choose the scale of as to be that of the momentum transfer at 

^ The most recent compilation of the Particle Data Group [18] suggests that the actual uncertainty in rric 
may be a factor of two smaller than the uncertainty that wc use here. However, since it is not clear that 
the systematic errors are well understood in the various determinations that enter into that compilation, 
we make a conservative choice of error bars. 
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the virtual-photon-charm-quark vertex, namely, mj/^. In the case of rjc — > 77, we choose 
the scale of as to be that of the momentum transfer at either of the photon-charm-quark 
vertices, namely, In order to take into account uncertainties in the scale and omitted 

corrections to the decay rates of next-to-next-to-leading order (NNLO) in a^, we attach an 
uncertainty to whose relative size is ctg. Then, we have 

as{mj/^) = 0.25 ± 0.06, (48a) 
as{m,jj2) = 0.35 ± 0.12. (48b) 

We choose the scales for the running QED coupling a to be the same as those for a^: 

Oi{mj/^) = (49a) 
(^irnj2) = (49b) 

where we ignore the uncertainties in a. 

In the case of the rjc matrix elements, we actually make use of two methods of computa- 
tion. One method is to compute the rjc matrix elements from r[r]c — > 77] , as we have outlined 
above. A second method is to equate the rjc matrix elements to the J/ ip matrix elements that 
we determine from T[J/ip e+e^]. Owing to the approximate heavy-quark spin symmetry 
of NRQCD [1], this equality is valid up to corrections of relative order v^. By combining 
these two methods of determining the rjc matrix elements, we can reduce the uncertainties. 
This approach is useful because the experimental result for r[77c — > 77] [Eq. (44b)] has a 
relative uncertainty that is comparable to the corrections to the spin-symmetry relation, 
which are of order ~ 30%. In principle, we could apply a similar approach to the J/ip 
matrix elements, but we would not gain a significant reduction in the uncertainties because 
the relative uncertainty in the experimental result for r[J/ijj e^e~] [Eq. (44a)] is small 
compared to the corrections to the spin-symmetry relation. In averaging the two sets of 
77c matrix elements, we must take into account the fact that many of the uncertainties are 
correlated between the two sets. We describe the procedure that we use for doing this in 
detail in the next section. 



^ We compute and a at each scale by making use of the code GLOBAL ANALYSIS OF PARTICLE 
PROPERTIES (GAPP) [23]. 
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B. Sources of uncertainties 



Let us now list the various uncertainties that enter into the calculations of the matrix 
elements. There is a theoretical uncertainty in the value of {q^) h that arises from the fact 
that the leading-potential approximation is accurate only up to corrections of relative order 
v^. For the computation that is based on T[J/ip —>■ e'^e~], we denote this uncertainty by 
A(qr^)j/,0, and for the computation that is based on T[ric — * 77], we denote this uncertainty 
by A{q'^)n^. We take these uncertainties to be v"^ ~ 30% times the central values. The 
uncertainties that arise from the scale uncertainties in ag and from neglecting NNLO cor- 
rections to the J/t/j and rjc electromagnetic widths are denoted by ANNLOj/^, ANNLO^^, 
respectively. As we have explained above, we parametrize these uncertainties as uncertain- 
ties in as [Eq. (48)]. However, we take ANNLOj/^ and ANNLO^^, to be uncorrelated. 
There are also uncertainties that are associated with the charm-quark mass rric [Eq. (45)], 
the string tension a [Eq. (46)], and the uncertainties in the experimental measurements of 
T[J/i/j e+e"] and r[rjc 77] [Eqs. (44a) and (44b)]. We denote these uncertainties by 
Amc, A(7, AFj/^, and AF^^, respectively. When we combine the values of the rjc matrix 
elements that we obtain from r[r)c — > 77] with those that we obtain from r[J/il) — > e'^e"] 
by invoking the heavy-quark spin symmetry, there is an uncertainty from corrections to the 
spin symmetry, which applies to the latter set of matrix elements. We take it to be f ^ ^ 30% 
times the values of that set of matrix elements. Since {q'^)j/^ already has an uncertainty 
A{q'^)j/^ of order v^, we apply this additional order-v^ uncertainty only to {Oi)j/^. We 
denote it by Ai>^. 

In making these uncertainty estimates, we have assumed that the standard NRQCD 
power-counting (velocity-scaling) rules [1] hold. Various alternatives to the NRQCD power- 
counting rules have been suggested [4, 24, 25]. Application of these alternative rules would 
affect our estimate of the correction to the heavy-quark spin symmetry, Af^, and our es- 
timates of the corrections to the static potential, A(g^)j/^ and A(g^)^^. In the standard 
NRQCD power-counting rules, Ai>^ is of relative order v'^. In the strong-coupling regime of 
Refs. [4, 24] and in the power-counting rules of Ref. [25], Av^ is of relative order Aqcd/"^c, 
which is actually smaller numerically than v^. The leading correction to the static potential 
is denoted by V^^^ /m (Ref. [4]). In the standard NRQCD power-counting rules, V'^^^ /m is 
suppressed as v"^ relative to the static potential. In the strong- coupling regime of Refs. [4, 24] 
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and in the power-counting rules of Ref. [25], V^^^m is of the same order as the static po- 
tential [24, 26]. In the lattice calculation of Ref. [27], V'^^ym corrects the string tension 
by about 17%, which is numerically smaller than v^. Other lattice calculations [28-30] also 
suggest that terms of higher order in the standard NRQCD power counting are suppressed at 
least as much as would be expected from the standard power counting. Therefore, we believe 
that the standard NRQCD power-counting rules give an upper bound on the uncertainties, 
and we use them for our uncertainty estimates. One could implement the alternative power- 
counting rules by equating Av'^ to Aqcb /ttIc times the central value and by by equating 
^(9^) J/V A(q^)^^ to 100% of the central value. 



C. Numerical results 



1. Computations using V[J/tp — > e+e ] and T[r]c — >■ 77] 

The results of our computations of matrix elements from T[J e+e^] and T[ric 77] 
are shown in Tables I and II, respectively. In each table, in the first row below the headings, 
we give the central values for the potential- model parameter A, the matrix element (Oi)if, 
and the ratio {q^)H- Subsequent rows contain the values for A, the matrix element, and the 
ratio that result from shifting each uncertain quantity in the calculation by plus or minus 
its uncertainty. We put a superscript 77 on the matrix element and the ratio for the rjc 
that are shown in Table II, in order to specify that these numbers are the result of a fit to 
^Vc 77]- 

The matrix elements and ratios of matrix elements, along with their uncertainties, are as 
follows: 

{0,)j/^ = 0.440lH?^lHJ^lH^^lH^^iHii GeV^ = 0.4401^1 GeV^ (50a) 

- 0.44ll°:}il°:°°^+°;°-l°:°^^+S;:i GeV^ = 0.44112;}!^ GeV^ (50b) 



//-OA 77 _ n 4Q4+0.042+0.083+0.012+0.083+0.112 q y3 _ r, 404+O.I69 Q yS /r-i N 

\'-^l/?7c —^■^'^^-0.040-0.069-0.012-0.066-0.118 '^®^ — ^■^>54_o^58 VjC V , {013,) 

/„2\77 _ r, ..0+0.133+0.024+0.038+0.023+0.041 r< y2 _ 440+O.I48 q y2 /r-,1 X 

\y /ric ~ '-'■^^'J-0. 133-0.028-0.037-0.028-0.038 — '-'■^^'^-0.149 ■ 

In the first equalities in Eqs. (50) and (51), the uncertainties are presented in the same order 
as in Tables I and II. In the last equalities in each of these equations, we have added the 
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uncertainties in quadrature. However it must be kept in mind for many applications that 
the individual uncertainties are correlated between the matrix elements. The correlations 
can be determined from the tabulations in Tables I and II. 

From Eqs. (10), (50), and (51) and the uncertainty in nic in Eq. (45), it can be deduced 
that 

{v')j/^ = 0.225^°:^°^, (52a) 
{vX = 0.22612;^^^ (52b) 

The central values of these results are somewhat smaller than an estimate, based on the 
NRQCD velocity-scaling rules [1], that (v^) should be equal approximately to v"^ ~ 0.3. 
However, they are consistent with being of order v^. 

We can see the effect of resummation by repeating our analysis, but keeping only the 
order-^;^ corrections in the formulas for the decay rates in Eqs. (22) and (31). The results 
are that the central values are shifted to (Ci) j/v- = 0.446269 GeV^ {q^)j/^ = 0.438520 GeV^, 
(0^)77 = 0.459867 GcV^ and {q^ = 0.433879 GeV^. Hence, the effects of the resum- 
mation on these quantities are —1.4%, -|-0.5%, —5.7%, and -|-2.1%, respectively. The small 
effects from resummation suggest that the v expansion of NRQCD converges well for the 
widths rfJ/'^ — > e+e~] and r[rjc — > 77]. 

2. Average values of rjc matrix elements 

Because some of the uncertainties in Tables I and II are correlated, we must take care in 
combining the results in these tables to obtain average values for the rjc matrix element and 
the Tjc ratio of matrix elements. First, we construct a two-by-two covariance matrix for the 
quantities and {Oi)'^'J from the deviations from the central values that correspond 

to the uncertainties listed in Tables I and II. Then, we use the inverse of the covariance 
matrix to construct for the deviation of the average value of (Oi),,^ from the two input 
values. We fix the average value of (Oi)?;^ by minimizing this with respect to it. The 
minimum value of is 8.9 x 10"''. This small value of refiects the fact that 
and are much closer in value than one would expect from the velocity scaling rules 

of NRQCD. Once we have obtained the average value of (Ci)r)c) we use it as an input to the 
potential model to compute the average value of {q'^)r]c- We carry out this computation of 
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the average values of (Ci)/,^ ^^'^ {q'^)vc values of the input parameters that correspond 
to each of the uncertainties that we have described. (The effect of the uncertainty A(qf^)^^ 
on the average value of (Ci),,^ has already been taken into account through the inputs to 
that average. We obtain the effect of A(q^)j^^ on the average value of {q'^)r^^ by varying the 
central value of the average value of (g^)r)c by ~ 30%.) The average values of (Oi),,^ and 
{Q^)vc that result from these computations are shown in Table III. The first row after the 
headings in Table III gives the central values of the averages of the rjc matrix element {Oi)r,^ 
and ratio of rjc matrix elements {q^)r]c- Subsequent rows show the effects of the various 
uncertainties on the average values. The central values and uncertainties in Table III can 
be summarized as follows: 

//n \ — n 4Q7+0. 024+0.033+0. 007+0.036+0. 006+0.073+0. 037+0.050 p„v3 — A'^T+O-Hl PoV^ 
\*^l/??c — -0.023-0.025-0.007-0.029-0.006-0.073-0.029-0.053 "'-'•^"^'-0.105 ' {00<X) 

/«2\ _ n 44O+0.132+0.009+0.040+0.010+0.002+0.026+0.010+0.018 /^o^/S _n/14O+0-143 r'^-\/2 /t;QK^ 
\y /ric — *-'-^^^-0.132-0.011-0.039-0.012-0.002-0.025-0.013-0.017 " *-'-^^^-0.143 ■ \'->'i'->) 

In the first equalities in Eq. (53), the uncertainties are presented in the same order as in 
Table III. In the last equalities, we have added the uncertainties in quadrature. As wc have 
mentioned, it must be kept in mind for many applications that the individual uncertainties 
are correlated between the matrix elements. The correlations can be determined from the 
tabulations in Table III. 

The correlated errors can also be expressed conveniently in terms of a correlation ma- 
trix. We construct a (symmetric) correlation matrix whose rows and columns correspond 
to {Oi)j/^, {q^)j/.,p, {^i)r)ci a-iid (9^))7c) respectively, taking the deviations from the central 
values from Table I and Table III. The result is 

/ 3.71 X 10-^ 1.64 X 10-^ 1.94 x 10"^ 8.18 x 10"^ 
1.64 X 10-^ 1.96 X 10-2 2.84 x 10"=* 1.93 x 10"^ 

1.94 x 10-3 2.84 x 10^3 hq ^ _3 7^ ^ ^q-a 

\ 8.18 X 10"^ 1.93 X 10-2 -3.71 x 10"^ 2.04 x 10"^ 

where the quantity in i-th row and j-th column is expressed in units of GeV"*'*'"'^ , with 
^1 = ''^3 = 3 and 712 = ^4 = 2. In charmonium decay and production processes, the 
NRQCD short- distance coefficients typically depend on rric- Hence, there may be correlations 
between the matrix elements and short-distance coefficients with respect to the uncertainty 
in rric. Therefore, we also give the correlation matrix for (Oi) j/^, {q^)j/ii;, (C'i)r?c) {Q^)vc^ 



(54) 
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(55) 



respectively, in which we omit the uncertainties that arise from mc- 

( 3.63 X 10-3 1.93 X 10-^ 2.21 x lO^^ 7.27 x lO^A 

1.93 x 10"^ 1.96 X 10-2 2.75 x lO'^ 1.94 x 10"^ 

2.21 X 10-^ 2.75 X 10-3 1.08 x 10"^ -8.86 x 10"^ 

\ 7.27 X 10-^ 1.94 X 10-2 -8.86 x 10-^ 2.03 x 10-^/ 

where the dimensions of the elements of C2 are the same as those of the corresponding 
elements of C\ in Eq. (54). We note that both correlation matrices C\ and C2 contain large 
off-diagonal elements that correspond to a correlation between the uncertainty in {q^)j/^ 
and the uncertainty in {q^)^^. Most of this correlation arises from the uncertainty in the 
string tension a. 



TABLE I: The potential-model parameter A, the NRQCD matrix element and the ratio 

{q^)j/^, as obtained from r[J/V' e^e~]. The first row below the headings contains the central 
values for A, the matrix element, and the ratio. Subsequent rows contain the maximum and 
minimum values for these quantities that are obtained by varying them with respect to each 
uncertainty. 



Case 


A 


(GeV3) 


{q^)j/^ (GeV2) 


central 


1.243 


0.440 


0.441 




1.256 


0.450 


0.573 




1.230 


0.430 


0.308 


+Amc 


1.233 


0.433 


0.443 


— Attt-c 


1.258 


0.451 


0.437 


+Aa 


1.191 


0.443 


0.482 


-Aa 


1.297 


0.437 


0.400 


+ANNLOJ/V, 


1.325 


0.504 


0.419 


-ANNLOj/^ 


1.166 


0.387 


0.459 




1.258 


0.451 


0.437 




1.228 


0.429 


0.444 
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TABLE II: The potential-model parameter A, the NRQCD matrix element {Oi)rf^, and the ratio 
{q'^)rjc, as obtained from T[r]c 77]. The first row below the headings contains the central values 
for A, the matrix element, and the ratio. Subsequent rows contain the maximum and minimum 
values for these quantities that are obtained by varying them with respect to each uncertainty. 



Case 


A 


(Oi);5:(GeV3) 


{q^nl (GcV^) 


central 


1.234 


0.434 


0.443 




1.291 


0.476 


0.576 




1.175 


0.393 


0.310 


+Amc 


1.340 


0.517 


0.415 




1.129 


0.364 


0.467 


+Aa 


1.195 


0.446 


0.481 


-Act 


1.276 


0.422 


0.406 


+ANNLO^, 


1.340 


0.517 


0.415 


-ANNLO^^ 


1.134 


0.368 


0.466 




1.374 


0.546 


0.405 


-AT,, 


1.041 


0.315 


0.484 



VI. COMPARISONS WITH PREVIOUS CALCULATIONS 

Our results for the matrix elements can be compared with those in Ref. [10]. In that 
paper, the values (Oi) j/^^ = 0.335 ± 0.024 GeV=^ and {Oi)^^ = 0.297 ± 0.032 GeV=^ are 
given. In the case of {Oi) j/,^, our result is 31% larger than that in Ref. [10]. Approximately 
6% of that change is the result of the change in the experimental value of V[J/il) — e^e~] 
from 5.26±0.37 keV [31] to 5.55±0.14±0.02 keV [18]. An implicit relativistic correction of 
about 22% arises from the use of mj/^ in Eq. (22), rather than 2mc. The use of a{rnji^) — 
1/132.6 [Eq. (49a)], rather than a — 1/137, decreases {Oi)j/^ by approximately 6%. The 
remaining change of about 9% is the result of including the explicit relativistic corrections 
in Eq. (22). In the case of our result in Eq. (51) is 46% larger than the value 

(C»i)]^L ^ Q ± 0.032 GeV^ that is given in Ref. [10]. In this case, there is a decrease in 
the value of (Oi)JJJ of 4%, owing to the change in the experimental value of T[ric 77] from 
7.5 ± 0.8 keV [31] to 7.2 ± 0.7 ± 2.0 keV [18]. The use of a{mr,j2) = 1/133.6 [Eq. (49b)], 
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TABLE III: Average values of the NRQCD matrix element (C'i)r?c the ratio (g^)r/c- The method 
of averaging is described in the text. The first row below the headings contains the central values 
for the matrix element and the ratio. Subsequent rows contain the maximum and minimum values 
for these quantities that are obtained by varying them with respect to each uncertainty. 



Case 


(Oi),, (GeV3) 


(GeV2) 


central 


0.437 


0.442 




0.461 


0.574 




0.414 


0.309 


+Amc 


0.470 


0.430 


—Arric 


0.413 


0.450 


+Aa 


0.444 


0.482 


-Aa 


0.431 


0.403 


+ANNLOJ/V, 


0.473 


0.429 


-ANNLOj/^ 


0.408 


0.452 


+Arj/^ 


0.443 


0.440 


-AFj/^ 


0.431 


0.444 


+Av^ 


0.511 


0.417 


-Av^ 


0.364 


0.467 


+ANNLO^, 


0.474 


0.429 


-ANNLO^^ 


0.408 


0.452 


+Ar^, 


0.487 


0.425 




0.385 


O.IGO 



rather than a — 1/137, decreases (Oi)]^ by approximately 5%. The use of as{mrj^/2) — 0.35 
[Eq. (48b)], rather than as{mj/^) = 0.25 [Eq. (48a)], which is used for the process rjc 77 
in Ref. [10], enhances the matrix element by approximately 14%. The remaining change of 
about 41% arises from the relativistic corrections in Eq. (31). 

In Ref. [10], the values {q'^)j/^ = 0.43 GeV^ and {q^)'Ji] = 0.25 GeV^ were obtained by 
making use of the Gremm-Kapustin [3] relation for the physical quarkonium mass and rric- 
While these results are not far from those in Eqs. (50) and (51), the uncertainties given in 
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Ref. [10] are on the order 100%, owing to the uncertainty in rric- In our calculation, we have 
been able to reduce the uncertainties significantly by making use of the Gremm-Kapustin 
relation (11) for the binding energy in the potential model. This leads to much smaller 
uncertainties than the use of the Gremm-Kapustin relation for the physical quarkonium 
mass and rric because we compute the binding energy directly in the potential model, instead 
of expressing it as a difference between niH and 2mc- 

In Ref. [2], the result {q'^)j/^ — 0.50 ± 0.09 ± 0.15 GeV^ was obtained from a potential- 
model calculation, which also made use of the Gornell potential. That result agrees, within 
errors, with the result in Eq. (50). In Ref. [2], the value of the matrix element {Oi)j/-^ was 
taken from Ref. [10], in which the relativistic correction to T[J/iIj e^e~] was not taken 
into account. The inclusion of that correction in the present work, along with a more precise 
determination of the potential- model parameter A, accounts for the difference in the value 
of {q'^)j/^ between Ref. [2] and Eq. (50). 

We can also compare our results with those in Ref. [32]. In that work, the following 
values are reported: {Oi)^J^ = 0.573 GeV^ (Oi)^J^ = 0.432 GeV^ and {Vi)^J^/ml = 
{Pi)^J^/ml = 0.0514 GeV^ These values were obtained by comparing the theoretical 
formulas for T[J/iJj — > e^e^], T[r]c ^ 77] i and T[J/iIj light hadrons] with the experimental 
results and by assuming that {Vi)^^'^ — {Vi)^J^. In the case of r[J/'0 light hadrons], 
processes involving an intermediate virtual photon were excluded in both the theoretical 
formula and the experimental rate. The theoretical expressions that were used in Ref. [32] 
to obtain these results contain the QGD corrections of relative order ag and the relativistic 
corrections of order v^. Taking rric = 1-5 GeV, which is the value that is used in Ref. [32], we 
find that the resuhs given in Ref. [32] yield {q^)^J^ = 0.202 GeV^ and {q^)^^^ = 0.268 GeV^. 
These values are considerably below the values in Eqs. (50) and (51) and considerably below 
the expectations from the velocity-scaling rules of NRQGD. The small values of (q^) and 
(q^)^J*-' are traceable to the use of the theoretical expression for r[J/j/j — > light hadrons]. 
In that expression, the coefficient of the contribution that is proportional to {q'^) j/^/ml is 
about —5.32 relative to the leading contribution. Because of this large negative coefficient, 
the quantity {q'^)j/^ must be much less than the values that wc obtain in order for the decay 
width to be positive. We regard this as an indication that the v expansion is not reliable for 
the rate r[J/ip light hadrons]. It is possible that the resummation methods that we have 
used in the present work could be used to tame the I'-expansion for r[J/'0 — > light hadrons]. 
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The value of {Oi)^J^ is about 30% larger than the value in Eq. (50) while the value of 
(Ci)^J*-' is about 1% smaller than the value in Eq. (51). Some of this difference is accounted 
for by the smaller values of (q^) j/v> and (9^),,^ in Ref. [32]. Ref. [32] also makes use of shghtly 
different values of rric (1.5 GeV) and (0.26) than those employed in the present work. A 
further difference is that the expressions for T[J/ijj — >■ e~^e~] and T[f]c — >■ 77] in Ref. [32] are 
expanded to first order in as and v^, rather than expressed as exact squares of amplitudes, 
as in Eqs. (22) and (31). 

Finally, there are quenched lattice computations [28] of the ground-state S-wme char- 
monium matrix elements that yield (Ci) j/^-^^ = 0.3312 ± 0.0006 ± 0.0030t[J:^5|:^ GeV^ and 
io^) j/'tl!-fic — 0.07-0.82 GcV^. In (Oi) j/,^-.,^^, the first error bar is from lattice statistics, the 
second error bar is from lattice systematics, and the third error bar is from the uncertainty 
in the one-loop perturbative computation that relates the lattice-regulated matrix elements 
to the continuum MS matrix elements. The lattice computations do not distinguish between 
the J/ip state and the rjc state. The lattice results are in agreement with our results, within 
uncertainties, but the lattice uncertainties are much larger than ours. These large uncer- 
tainties arise from the uncertainty in the perturbative conversion from lattice to continuum 
MS matrix elements. 

VII. SUMMARY 

For many jS'-wave heavy-quarkonium decay and production processes, the color-singlet 
5'-wave NRQCD matrix elements of leading order in v enter into the dominant theoretical 
contribution. The first relativistic corrections to these processes involve the matrix elements 

of relative order v^. 

We have computed the color-singlet 5'-wave NRQCD matrix elements of leading order 
and next-to- leading order in v"^ for the J/ip and the r}c- For each of these quarkonium 
states, we have determined the values of these matrix elements by comparing the theoretical 
expressions for the electromagnetic decay rates (r[J/'0 — > e+e"] or r[77c — > 77]) with the 
experimental measurements and by using a potential model to compute the matrix elements 
of relative order v"^. If the static, spin-independent QQ potential were known exactly, then 
the potential-model calculation would be accurate up to corrections of relative order v^. We 
made use of the Cornell potential and fixed its parameters by using as inputs the lattice 
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measurements of the string tension, the J/ip-ip{2S) mass sphtting, and the quarkonium wave 
function at the origin, which corresponds to the NRQCD matrix element of leading order in 
V. Because the potential-model calculation of the order-v^ matrix element depends on the 
leading-order NRQCD matrix element and the decay widths depend on both of these matrix 
elements, we obtained the matrix elements for the J/ip and the rjc by solving, in each case, 
two coupled nonlinear equations. 

In the theoretical expressions for the electromagnetic decay widths, we made use of the 
generahzed Gremm-Kapustin relation (11) (Ref. [2]) to resum a class of relativistic correc- 
tions. This resummation includes all of the relativistic corrections that are contained in the 
leading-potential approximation to the quarkonium QQ color-singlet wave function, up to 
the ultraviolet cutoff of the NRQCD matrix elements. 

There are many sources of uncertainties in our calculation. Some of these are correlated 
among the matrix elements. Therefore, we reported the variations of the matrix elements 
with respect to each source of uncertainty. 

The experimental measurement of the width r[r]c — > 77] has relatively large uncertainties, 
which translate into large uncertainties in the r)c matrix elements. Owing to the heavy-quark 
spin symmetry [1], the J/tjj and rjc matrix elements arc equal, up to corrections of relative 
order v"^. Therefore, we were able to reduce the uncertainties in the rjc matrix elements by 
averaging the values that we obtained from r[r)c — > 77] plus the potential model with the 
values that we obtained from r[J/ip — > e'^e"] plus the potential model. In performing this 
average, we took into account the additional uncertainty of relative order v'^ that arises from 
equating r/c matrix elements to J/i(j matrix elements. 

Our principal results are given in Tables I, 11, and 111 and are summarized in Eqs. (50), 
(51), and (53) for the matrix elements that were determined from r[J/i/j e"'"e~], the matrix 
elements that were determined from r[r)c — > 77], and the average of the two, respectively. We 
consider the results in Table III and Eq. (53) to be our best values for the 77c matrix elements. 
In applying these results to calculations of quarkonium decay and production rates, it should 
be kept in mind that the uncertainties are highly correlated between matrix elements and 
that there arc correlations between matrix elements and short-distance coefficients with 
respect to the uncertainties in rric- Therefore, it may be necessary to use all of the information 
that is contained in Tables I, II, and III, rather than to rely on the summaries in Eqs. (50), 
(51), and (53). 
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Our results in Tables I and II and Eqs. (50) and (51) conform to the expectation, from 
the heavy-quark spin symmetry, that the J/%1) and r]c matrix elements are equal, up to 
corrections of relative order ^ 30%. In fact, the leading-order matrix elements differ by 
about 1.5%, while {q^)j/^ and {(f'T^^ differ by only about 0.5%. The velocity-scaling rules 
of NRQCD [1] state that the quantities {v'^)h = {Q'^)H/n^l should be of order v'^ ?a 0.3. 
From Eq. (52), it can be seen that our results satisfy this expectation, although they are 
somewhat smaller than the nominal value of v'^. 

As we have discussed in Sec. VC 1, the effects from resummation on our results are small, 
ranging from —5.7% for (Oi)^J to 2.1% for {q'^)^']- The small effects from resummation 
suggest that the v expansion of NRQCD converges well for the widths T[J/ip — >• e+e"] and 

rfe ^ 77]- 

Our results for and {Oi)^'^ are considerably larger than those in Ref. [10], pri- 

marily because we have included relativistic corrections to the electromagnetic decay rates 
in the present work. The changes in the values of these matrix elements would significantly 
increase the rate for the process e+e" — > J/t/j + rjc that is calculated in Ref. [10]. 

Our result for {q'^)j/^ agrees, within uncertainties, with that in Ref. [2], but is slightly 
smaller. Most of this difference arises from the fact that, in Ref. [2], the value of 
was taken from Ref. [10]. 

In Ref. [32], a much smaller value for {q^)j/^ was reported. {{q^)'q^ was assumed to be 
equal to {q^)j/iij in this work.) The smallness of {q^)j/ijj in Ref. [32] can be traced to the 
use of the width r[J/'0 — > light hadrons] to constrain the matrix elements. The theoretical 
expression for that width contains large order-v^ corrections that, in our opinion, make 
the reliability of the expression suspect. It is possible that the resummation technique 
that we have employed in this paper could be used to bring the v expansion for r[J/'0 
light hadrons] under control. 

Our results for {Oi)j/^ and {q^)j/.^ are in agreement with those from lattice calculations 
[28] , although the lattice uncertainties are much larger than ours. 

We beheve that the values that we have obtained for the J/ip and rjc color-singlet NRQCD 
matrix elements are the most precise ones that are available to date. The new values for the 
matrix elements of leading order in v should have a significant impact on the calculations of 
a number of charmonium decay and production processes [33] . For quite a few charmonium 
processes, it is clear that relativistic corrections are important. Within the framework 
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of NRQCD, the matrix elements of order are essential ingredients in calculating those 
corrections. We have also attempted to quantify all of the significant theoretical uncertainties 
in our determination of the J/ip and r^c color-singlet matrix elements. Our treatment of 
uncertainties could provide the basis for more reliable estimates of theoretical uncertainties 
in future calculations of charmonium decay and production rates. 
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